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1. Introduction 

When jets are produced at colliders, they can be initiated either by a gluon or by a quark. Be- 
cause of the fact that these two partons carry different color charges and spin, it is expected that the 
corresponding jets show different properties: typically a gluon jet is braoder and contains a larger 
amount of hadrons. Jets with different parent partons can also be distinguished and investigated 
looking for the jet charge distribution as discussed in detail in Ref.[[I]] with important applications 
for the LHC. The interactions between quarks and gluons are described by quantum chromodynam- 
ics (QCD) and to understand the difference between a quark and a gluon jet is one of the strongest 
test of this theory. 

The typical way to depict the production of a jet from a parton is the following: an initial quark 
(or gluon) start radiating gluons, which in turn can radiate further gluons or split into secondary 
quark-antiquark pairs. In this way the virtuality of the parent parton decreases in a so called parton 
showering process. Finally when the virtuality falls below a certain cutoff the cascade stops and 
the final state partons hadronize into color neutral hadrons, a process usually described by phe- 
nomenological models. This happens because the production of hadrons is a typical process where 
non-perturbative phenomena are involved. However, for particular observables, this problem can 
be avoided. In particular the "counting" of hadrons in a jet which is initiated at a certain scale Q 2 
belongs to that class of observables and in this case one can assume with quite high accuracy the 
hypothesis of Local Parton-Hadron Duality (LPHD) which simply states that parton distributions 
are just renormalized in the hadronization process without changing their shape This would in 
principle allow perturbative QCD at fixed order to make predictions without the need to consider 
phenomenological models of hadronization. Nevertheless these are processes dominated by soft 
gluon emissions and it is a well known fact that in such kinematic regions of the phase space fixed 
order perturbation theory fails and resummation is needed (see e.g. Ref.[|3|]) and this is the main 
topic discussed in this talk. 

The gluon and quark multiplicities (^(Q 2 )}^ and {nh(Q 2 )) s represent the avarage number 
of hadrons in a jet initiated by a gluon or a quark respectively at scale Q 2 . In the past analytic 
predictions have been achieved solving the equations for the generating functionals in the modified 
leading logarithmic approximation (MLLA) in Ref.[Q] up to the so called N 3 LOr in the expansion 
parameter y/c^ i.e. a s . However for the ratio r = (nh) g / (n/t)j the theoretical prediction is about 
10% higher than the data at the scale of the Z° vector boson and the difference with the data 
becomes even larger at lower scales even if the convergence of the perturbative series behaves very 
well. An alternative approach was given in Ref. [Q] where equations for the derivative of the ratio 
of the multiplicities are obtained in the MLLA within the framework of the colour dipole model. 
There a constant of integration which is supposed to encode non-perturbative contributions is fixed 
by the data. A constant offset to the quark and gluon multiplicities has also been introduced in [§]. 

Very recently a new formalism has been proposed and developed in Refs.[j7|, g], where the 
problem of the apparent good convergence of the perturbative series is solved and where any ad 
hoc offset is needed once one includes the effects coming from the full mixing between quarks 
and gluon evolutions. The result looks like a generalization of the result obtained in Ref.[Q]. In 
the new approach the non-perturbative physics to the gluon-quark multiplicities is parametrized 
simply in the initial conditions of the evolution equations. Due to the good agreement with the 
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data obtained in Ref.[j8]], at the time of the contribution to this conference we have done a study 
where we have extended our analysis adding as a fit parameter also the strong coupling constant. 
Thanks to the new outstanding results in small-x timelike resummation obtained in [§] in the MS 
scheme, we were able to compute exact next-to-next-to-leading-logarithm (NNLL) contributions 
to the evolution of the multiplicity with with approximated NNNLO normalization factors in the 
yfa~ s expansion. Previously published results where available up to the NLL accuracy but in a 
massive gluon scheme which in general is unfortunately not suitable to combine resummation with 
fixed higher order corrections which are naturally given in the MS scheme. See Refs.fliOl, [TT| ] for a 
general discussion about the scheme choice and scheme dependence in this context. 

In the following after reviewing the history which lead to the latest important improvements in 
"timelike" QCD of Refs.lO, 0] and sketching our formalism, we present, as anticipated, the result 



of a slightly different global fit than in [g] where also the strong coupling contstant has been added 
among the free parameters confirming the good quality of the fit. 



2. Exiting times for "timelike" QCD 



The important and fundamental results obtained in Refs.[ 12, H] represent the beginning of the 



happy end of a quite long story: In the 1972 it was realized in Ref.[|13[] that at the lowest order the 
"timelike" splitting functions (occuring in Semi-Inclusive electron positron Annihilation (SIA)) are 
the same as their "spacelike" counterparts (occuring in Deep Inelastic Scattering (DIS)). This prop- 
erty goes under the name of Gribov-Lipatov relation. Then in the 1980 it has been shown by Curci, 



Furmanski and Petronzio [14] that the Gribov-Lipatov relation is violated. In that period several 
groups obtained the NLO contributions to the 'timelike' splitting functions [|14[ ]B[ |T^, [17|, |lj| [Pj]. 
In 2004 the NNLO "spacelike" splitting functions were published by Moch, Vermaseren and Vogt 



in Ref . []20j, [21| ]. In this last case the calculation could be performed via forward scattering ampli- 
tudes, a fact that has allowed a direct calculation in terms of Feynman diagrams. This is however 
not the case for the "timelike" splitting functions and different techniques have been investigated. 
Two of them have been the most successfull. The first one is due to Dokshitzer, Marchesini and 
Salam [|2|] who developed a formalism trying to rescue (at least at the formal level) the Gribov- 
Lipatov relation at higher orders. The approach of Ref. p2|] shaded light on many theoretical aspects 



revealing (using their words) "intrinsic beauty of the perturbative quark-gluon dynamics" [23p. The 
second one is based on the fact that the "timelike" and the "spacelike" splitting functions can be 
related by the analytic continuation of the scaling variable x — > l/x with x representing the fraction 
of the parton longitudinal momentum. In particular, even if beyond the LO this cannot be done 
directly from the splitting functions, it can be done for the corresponding physical evolution ker- 
nels and this has been shown in Refs. [|4], 25]. Formally, the evolution of a structure function 
is governed by a physical kernel K which depends on the coefficient function C and the splitting 
function P according to: 



Q- 



dF 



K®F 



0(*)£+c. 



P ®C 



-l 



where j5 is the QCD beta function for the running of the strong coupling a s 

,da s (Q 2 ) 



dQ 2 



(2.1) 



(2.2) 
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and where ® means the usual integral convolution with respect to x. Eq.(2.1) also tells us that 

F = C®D, (2.3) 

where D is what we call fragmentation function (the analog of the parton densities in the "space- 
like" case) and that 



Q 1 



dD 



P®D. 



(2.4) 



According to Eq.( 2. 1 ) and the discussion above it, we see that in principle all the "timelike" splitting 
functions can be obtained from the "timelike" physical kernels (obtained by the analytic continua- 
tion x — > 1 /x), once the "timelike" coefficient functions C are known. The "timelike" coefficients 
functions are known at NNLO 28. |2^, 3(J 31, 32], while the analytic continuation of the 



physical kernels from the NNLO "spacelike" ones has been obtained in Refs. []33|, |34[ ] and finally 
in Ref.[[l2|]. Here we mention only that the analytic continuation jc — > 1 /jc becomes subdle (see 
Ref.[25]) for the logarithmic terms singular in x — > 1 for which 



ln(l 



ln(l 



lnx + /7T, 



(2.5) 



and that additional constaints (see Ref . JT2| ] and referencies therein) from the momentum sum rules, 
the supersymmetric limit (Ca = Cp = n/) and from the relations found in Ref.[^2|] are needed to 
close the problem in a satisfactory way. 

Nevertheless the NNLO "timelike" splitting functions for the quark-singlet system obtained 



in Ref. [ |12| ] are not the end of the story. Indeed they present singularities in the threshold (x — > 
1) and in the large energy limit (x — > 0). Both kind of singularities are due to the radiation of 
soft gluons and they make perturbation theory to fail. As already mentioned, the difficulty with 
these singularities are overcome by resummation. For very recent developments in the large-x 



resummation see e.g. Refs. [[35], [36J]. Here we are mainly interested in the small-x resummation 
because it is directly related to the computation of hadronic multiplicities as already realized a long 
time ago [B]. With respect to the multiplicity studies, the basic equation is the one governing the 



evolution of fragmentation function D for the gluon-singlet system Eq.(2.4), which in Mellin space 
sees the convolutional product <g) turn into an ordinary matrix product: 



Q 



; d ( D s ((0,Q 2 ) 

d<2 2 \D g ((o,Q 2 ) 



P aq (co,a s ) P gq (co,a s 
Pq g ((Q,a s ) P gg (co,a s ) 



D s (co,Q 2 ) 
D g {(0,Q 2 ) 



(2.6) 



Here CO =N — 1 with N the standard Mellin moments with respect to x and where Py are splitting 
functions as introduced in Eq.(^7j]). The standard definition for the hadron multiplicity in terms of 



the fragmentation function is given by (see e.g. Ref. [37]) the integral over x of the fragmentation 
function which is just the first Mellin moment (ft) = 0): 



(n h (Q 2 ))a 



f 1 dxx a D a (x,Q 2 

Jo 



D a (co = 0,Q z 



(2.7) 



<u=0 



where a = s for a quark jet and a = g for a gluon jet. It is clear from the definition in Eq.(2.7) 
that eventual non-integrable singularities in x ~ in the inegrand could not occur for the definition 
to make sense. We remind that the singularities in x = are translated into singularities in co = 
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0. We have already mentioned that resummation which includes the singularities from all orders 
according to a certain logarithmic accuracy is the standard way to solve this problem. Indeed 
after resummation, the singular behavior in x ~ (or equivalently in ft) ~ 0) disappears. The 
generally better choice of the MS scheme for to perform resummation has been considered and 
solved only quite recently. Firstly the NLL accuracy has been obtained in Refs.[38, 3^] and finally 
the NNLL accuracy has been reached in Ref.fP]. Here thanks to the approach used by A. Vogt 



in Ref.jp 8[], the singularities of the splitting functions are iteratively extracted according to the all 
order factorization into the transition function Z(e) of the £ = poles in dimensional regularization 
(d = 4 — 2e). Indeed, according to the factorization theorem p0| , 14], we can rewrite Eq.2.3 in 
Mellin space as 

F = CD = C°(e)Z- l (£)Z(e)D , (2.8) 

where C° and D° are the "bare" coefficient function and fragmentation function. Z(e) is the tran- 
sition function containing only poles in e = that are factored out from C°. Hence, substituting 



Eq.(2.8) into Eq.(2.4) we get that the splitting functions can be directly related to the transition 
function in the following way: 



P = Q 2 



dZ 



7-\ 



o, \ 3Z 



(2.9) 



Now it is shown in Ref.[^] how one can solve this equation in Z obtaining at all orders the three 
highest order poles in e knowing the NNLO corrections to P and j8. Additionally knowing the 
higher order corrections to C which is pole free one obtains from the NNLO computations the all 
order structure of the three first highest singularities in e : 



C°(e)=C(e)Z(e). 



(2.10) 



The key point of Ref.jp 8[] is that they realized that for example for the case of the gluon the small 
ft) behavior of the bare coefficient function is 



C°(e) 



1 y Og 
ft) *-? s 2 "- 1 



n-1 



1 



■2(n-l)e/co 



(A ln + eB ln + e 2 C ln + ...)■ 



(2.11) 



Finally comparing Eq.(2. 11) with Eq.( 2.10| ) one obtaines sytems of equations for the coefficients 
A,B and C, which produce sequences up to arbitrary orders in a s of the three highest powers in 
1/ft) or equivalently (back to *-space) in hue. Then the highly non trivial part of this approach is the 
solution of these sequences that are obtained for the large logarithms and this has been successfully 



obtained in Refs.[38, 



3. Plus-minus component evolution for the gluon-singlet system 



In Ref.[]8|] it has been shown that, by use of the results obtained in Ref.[^], the scale dependence 
on Q 2 of the gluon and quark multiplicities defined by Eq.(2.7) and governed by Eq.(^) can be 
computed analytically in a closed form. The approach that we used is a generalization of the 
techniques used in Ref.[41] used to solve the evolution equation of the parton densities in the 
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"spacelike" case. Our result can be written in the following way 

{n h {Q 2 )) g = (n h (Q 2 Q )) g f + (Q 2 ,Q 2 ), 



{n h {Q l )) s = (n h (Q 2 )) 



2^ T + {Q\Ql 



r + (Q 2 ) 



+ 



(n h {Ql))s 



(nh(Q 2 )) g 
r + {Ql) 



?-(Q 2 ,Ql), (3.1) 



where <2o is an arbitrary reference scale (see below). The dependence on Q 2 is always inside the 
running coupling constant a s (Q 2 ). Indeed substituting the QCD values for the color factors and 
choosing rif = 5 in the formulae given in Ref.[||] we can write that at NNLL 



T-{Q\Ql) 



«,(6 2 ) 
a s (Q 2 ) 



exp < d2 



VocM 2 



+ d 3 Joc s (Q 



where 



di =0.38647, d 2 = 2.65187, ^ = -3.87674, d 4 = 0.97771. 



(3.2) 



(3.3) 



The function r + (Q 2 ) in Eq.dOl) has been computed in @, up to the third order (NNNLO) 
in the yfa~ s expansion and in QCD is given by: 



r + (0 2 ) = 2.25-0.61567 J a s {Q 2 ) - 2. 19156 a s (Q 2 ) + 0.24348 a, 3/2 (Q 2 ) + 0(a 



(3.4) 

again with rif = 5. The first term in Eq.(^4j) corresponds to the lowest order expectation for the 



ratio between the gluon and the quark multiplicities given by Ca/Cf- According to Eq.(3.1), our 
prediction for the multiplicity ratio r is 



r + {Q 2 ) 



l i r + (Q?) an h (Ql)) s r + (Ql) 



A T-{Q\Ql) 
' t + (Q\Ql) 



(3.5) 



From this expression we can see that our result is a generalization of the result in Eq.( |3.4[ ) due to 
the inclusion of the contribution proportional to f_ (Q 2 , <2g) in the denominator of Eq.(3.5). 



Now we show that our results in Eqs.(3.1,33|) do not depend on the reference scale Qq and 
doing this also the meaning of f±(Q 2 ,Ql) will become clear. We define the so called "plus" and 



"minus" component in terms of the singlet and the gluon components defined in Eq.([2.7|) in the 
following way: 



(n h (Q 2 )) + ee (n h (Q 2 )) g , 
(n h (Q 2 ))_ ee (n h (Q 2 )) s - 



MQ 2 )) g 
r + (Q 2 ) 



0-6) 
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With this definition we have that the system of Eqs.(3.1) is equivalent to the following one 



(n h (Q 2 ))± = {n h (Q 2 Q ))±f ± (Q 2 ,Ql). 



(3.7) 



By noting that according to Eqs.(3.2), 



t±{Q 2 ,Ql)t ± (QlQ\) 



(3.8) 



one can directly check that, after expressing {rih{Q 2 ) ))± in terms of Q\ using Eqs.(pV7|) and then 
replacing them again into Eqs.(3.7), one arrives at the same expression but with <2q replaced by 



Q 2 . This shows the independence of the result on the reference scale <2q which is arbitrary and 
that 7±(<2 2 ,<2o) just represent the renormalization group exponents in a basis where the splitting 



function matrix of Eq.(2.6) can be considered diagonal. To look for a basis where the splitting 
function matrix could be considered diagonal by a sufficient accuracy was the main motivation 
under the computation in Ref.[||]. 



4. Global fit to the avarage multiplicity data: determination of a s 



Our result Eqs.(3.1) depends on three quantities that should be fixed by the comparison with 
the experiment. They are {nh{Qfy) g , (nh{Qo)) s (representing the avarage multiplicity at the ref- 
erence scale <2q) and a reference value for the running coupling constant which is conventionally 
chosen at the scale of the mass of the Z° weak boson, a s (Mz). In Ref.[^], a s (Mz) has been fixed 
as an input parameter to be equal to the world avarage value Q]. The result that we obtained was 



(n h (Q 2 Q )) g = 24.02 ±0.36, (n h (Q 2 Q )) g 
oc s (M z ) = 0.118, input parameter, 



15.83±0.37, 90% C.L., 



(4.1) 



where the arbitrary scale has been chosen to be <2o = 50 GeV. We have checked that the % 2 remains 
unchanged changing the value of <2q as expected from the considerations done above. 

Due to the quite good agreement with the data obtained in Ref.[|8|], we tried to treat a s {Mz) as 
an additional fit parameter. The result that we obtained for such a fit (always for <2o = 50 GeV) is: 



(n h {Ql)) g 
a s (M z ) 



24.18±0.32, {n h (Ql)) g = 15.86 ±0.37, 90% C.L. 
0.124 ±0.005, 90% C.L., 



(4.2) 



which is in agreement with Eq.(|4.1|) and with the experimental values for the gluon and quark mul- 
tiplicities at 50 GeV. The % 2 per degree of freedom is improved in this fit from 3.71 (corresponding 
to Eqs.(4T)) down to 2.84 (corresponding to Eqs.(|4~2|)). The predicted value for a s (Mz) is higher 



than the world avarage value [44], however such a higher value is not new in hadronic final states 



of electron-positron annihilation [|45|]. However how much it would make sense to compare our 
result with af s is still under investigation. For the determination of the result in Eqs j4.2| we have 
used the same data selection we used in Ref.[||], which is also the one adopted by the DELPHI 
colaboration in Ref.[46]. In Fig.[j] we plot the result of our new fir according to Eq.(fO|) together 



with the uncertainties and the data, which are taken from the Tables in Ref.[47] 
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Figure 1: Gluon and quark multiplicities fits according to Eq.(4.2) compared to the data. 

5. Conclusions 

In this talk we have sketched the history and the various techniques culminating in the very 
important new results from the theoretical side (see Refs. [|12|, |9|] and referencies therein for the 
most recent developments), which have already made possible new important improvement in the 
descriprion of the data in hadronic final states processes in the "timelike" case (see Ref.|j8]]). Mo- 
tivited by the well description of the data obtained in Ref . [j|] for the avarage multiplicities in jets 
initiated by gluon and quarks, we have extended our previous analysis making the value of the 
strong coupling constant at the reference scale of the mass of the Z° boson a parameter for the 
fit. In this way we have still obtained numbers in agreement with observations for the value of the 
multiplicities at the reference scale <2o = 50 GeV and in addition we got a significant improvement 
in the % 2 per degree of freedom and, as a byproduct, a prediction for a s {Mz). 
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